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We study the boundary conditions at a fluid-solid interface using molecular dynamics simulations 
covering a broad range of fluid-solid interactions and fluid densities, and both simple and chain- 
molecule fluids. The slip length is shown to be independent of the type of flow, but rather is related 
to the fluid organization near the solid, as governed by the fluid-solid molecular interactions. 

PACS numbers: 51.10.+y, 34.10.+X, 92.20.Bk 



The principal theme of this paper is a study of the na- 
ture of the boundary conditions (BC) of fluid flow past a 
solid surface, a crucial ingredient in any continuum fluid 
mechanical calculation. The BC cannot be deduced from 
the continuum differential equations themselves, and it is 
often not easy to determine them experimentally. While 
the normal component of the fluid velocity must vanish 
at an impermeable wall for kinematic reasons, the paral- 
lel component, when extrapolated toward the wall, may 
match that of the wall only at some distance £ away from 
it. This phenomenon is known as slip and £ is the slip 
length 0. Since the pioneering work of Maxwell it 
has been recognized that the scale of ( for a simple dilute 
gas is set by the mean free path A of the fluid molecules, 
with an 0(1) proportionality constant for a thermaliz- 
ing wall. However, for a specularly reflecting wall, the 
proportionality constant could become large and lead to 
large slip. In the limit of low fluid density, p, A be- 
comes large suggesting that ( would be large as well. 
Furthermore, in this limit the continuum approximation 
need not hold, and it is hard to glean the nature of the 
BC without a detailed knowledge of the influence of the 
fluid-solid interaction ||. This is indeed the situation 
for micro-electro-mechanical systems 0] which operate 
in the Knudsen regime, in which A can be larger than 
the system size. Earlier molecular dynamics (MD) stud- 
ies have indicated substantial velocity slip for repulsive 
walls and deviations in hydrodynamic velocity profiles 
near the wall on lowering p 



In our MD simulations, we find that the flow profile in 
the middle of a channel does indeed correspond to that 
predicted by continuum theory, but we observe a range of 
behaviors near the walls. Our work provides a molecular 
basis for the large variability in the amount of slip ob- 
served experimentally Q . We find that £ is an excellent 
descriptor of the boundary conditions, independent of the 
channel width and the nature of the flow. Even at high 
densities, significant slip is induced on weakening the 



wall-fluid attraction. In the low p subcontinuum regime, 
a large £ is found in virtually all cases, except for a chain 
molecule liquid and a strongly attractive wall. These two 
distinct classes of behavior lead to predictions amenable 
to experimental test. First, in a pressure driven flow, the 
speed with which fluid is transported shows a maximum 
as p is varied for all the large slip situations and none 
when the slip length remains small. Second, this maxi- 
mum speed should scale linearly with the channel width 
for the large-slip case, but should scale quadratically (as 
in a dense fluid) for the small slip situation. Finally, our 
results indicate a strong link between the degree of slip at 
low p and the nature of fluid organization in the vicinity 
of the channel walls. In addition, we find analogous be- 
havior in calculations of thermal slip, which arises when 
the two bounding walls of a channel are held at different 
temperatures [jlO[ . 

In order to elucidate the BC as one lowers p from the 
dense liquid to the dilute gas regime, we consider a chan- 
nel geometry and adopt the parametrization of Thomp- 
son and Robbins ||. There are two parallel walls in 
the x — y plane and N fluid atoms between the walls. 
Periodic BC are imposed along the x and y directions. 
In the monatomic case, two fluid atoms separated by 
a distance r interact with the Lennard-Jones potential 
V LJ {r) = 4e [(^r 12 - (^r 6 ] , where a is the size of 
the repulsive core. As in Ref. the potential is trun- 
cated at 2.2cr (and shifted). The MD studies are per- 
formed at a temperature T — X.le/kB, which is just 
above the liquid-gas coexistence region for this model. 
The characteristic MD time unit tq = ^/ma 2 /e, where m 
is the atomic mass, is 0(1 ps) for typical fluids. The walls 
are constructed from two [001] planes of an fee lattice (96 
atoms at each wall, with lattice constant 0.85c), with 
wall atoms tethered to fixed lattice site by a harmonic 
spring with a large spring constant. For the narrowest 
channel studied here, the fluid atoms were confined to a 
volume 13.6cr x 5.1cr in the x — y plane and Lq — 12.75(7 



between the inner faces of the walls. The wall- fluid in- 
teractions are modelled by a distinct Lennard- Jones po- 
tential V wf {r) = 16e [(f)" 12 - A (£ )- 6 ] . The param- 
eter A determines the properties of the wall and varies 
between I and 0, corresponding to attractive and repul- 
sive walls, respectively. Our earlier studies [|llj of colli- 
sions between fluid molecules and this wall have quan- 
titatively confirmed that A=l corresponds to a purely 
thermal wall, and A=0 to a purely specular one. 

The chain molecules comprising the polymeric fluid are 
composed of n — fO Lennard- Jones atoms each. The 
consecutive atoms along the chain are tethered by the 
FENE potential |f], Vfene = -K,/2log[\ - (r/r ) 2 ], 
where re = 30e and tq = 1.5a, which together with Vlj 
keeps these atoms around r = a apart. This system is 
studied at T = 1.6e/ks, which is above the liquid-vapor 
critical point, and the molecular diameter is small com- 
pared to the channel width (the radius of gyration of 
a molecule does not exceed 1.7a). For either fluid, the 
equations of motion were integrated using a fifth-order 
predictor-corrector algorithm |f 3|] . The spatial averaging 
was carried out in slabs of width cr/4 parallel to the pr- 
axis. The high dilution data points were time averaged 
over 3 million To to reduce the noise. The results given 
here were obtained with a Langevin thermostat [fl3[ , but 
a Nose-Hoover thermostat [[13| was found to yield virtu- 
ally identical data. The T-dependent noise was applied 
in all directions during equilibration and only in the y- 
direction when data were collected. In the studies of 
thermal slip, the thermostat is applied only to the wall 
atoms, where we assign T = Lie/fee to one wall and 
T = 1.3e/kb to the other, and determine the variation of 
T across the channel. 

The density profiles for the thermal slip simulation of 
the monatomic fluid are shown in Figure 1, for the at- 
tractive wall case (A=l). One can identify several layers 
near the wall, but at the center of the channel the profile 
becomes smooth. The value of density in the center, p, 
is used in the remaining plots as a measure of the fluid 
density. (A layered density profile is also characteristic 
of both monatomic and polymeric fluids without an im- 
posed temperature gradient The density profile is 
asymmetric, as illustrated by the behavior of the density 
difference in the second layers between the cold and hot 
sides, Ap2 = p% c ~ P%hi shown in the inset of Figure 1. 

For both kinds of fluids, with and without a tempera- 
ture gradient, the density in the first layer is largely p- 
independent, as shown in the top two panels of Figure 
2 for the uniform-T case. There is, however, a notice- 
able difference in the behavior of the second layer. In 
the monatomic case, the density in the second layer un- 
dergoes a steady buildup from a negligible value, when 
p=0.004<7 -3 (corresponding to 7V=100, where 96 atoms 



coat the walls and the remaining 4 move essentially bal- 
listically within the channel), to a saturation value when 
p becomes large. In the polymeric case, on the other 
hand, the second layer is almost fully constructed even at 
p=0.002a~ 3 and then its density varies only in a minor 
way. Thus the effective boundaries for the monatomic 
fluid undergo rapid transformation as a function of p, in 
contrast to the behavior of the chain molecule fluid. 

The top and bottom panels of Figure 3 describe viscous 
and thermal slip, respectively, for the monatomic fluid. 
The velocity profile of the top panel was obtained for 
Couette flow. The central region has the expected hy- 
drodynamic behavior - a linear velocity profile with con- 
stant shear stress for Couette flow, and a parabolic veloc- 
ity profile with a linear shear for gravity-driven flow (not 
shown). While the temperature profiles in the thermal 
gradient case are remarkably similar to the velocity pro- 
files in Couette flow, the dependence of the thermal slip 
length (defined analogously to the velocity slip length) 
on p is non-monotonic in a way that reflects the behavior 
of Ap2- As seen in the middle panel of Figure 4 for A=l, 
the slip is a minimum at p as 0.2er~ 3 , just where Ap2 is 
maximal (see inset of Figure I). Generally, the profiles 
predicted by continuum theory are seen either to persist 
to a point close to the wall (no slip), or else to terminate 
sooner, with a rapid change in profile around the second 
layer. It is the latter situation which is associated with 
slip, and it arises at very low p in most cases. Recently, 
Travis and Gubbins Q] have performed MD studies of 
velocity profiles for A — j in gravity driven flows in a dif- 
ferent parameter range: liquid densities, much narrower 
wall separation 0(4ct), and with an acceleration an or- 
der of magnitude larger than used here. Presumably due 
to these differences, these authors obtained triple peaked 
velocity profiles with large contributions to the flow from 
the first layers. 

The precise definition of the slip length £ is given in terms 
of the position, zq, at which the extrapolated fluid ve- 
locity reaches that of the wall (non-zero in the Couette 
case) , as measured relative to the inner layer of the wall 
atoms. C is taken to be positive for z outside of the 
channel and negative otherwise. Figure 4 summarizes 
our results on viscous and thermal £ for monatomic and 
chain molecules. Viscous slip was studied in Couette and 
gravity driven flows and at various channel widths (two 
shown). We find a large variation in boundary conditions 
on changing the material parameters, in agreement with 
experimental observations ||. However, the value of £ 
is found to be transferable in that it neither depends on 
the kind of flow nor on the channel width. For a range 
of high p values for the monatomic fluid with wall in- 
teraction A=l, £ is roughly equal to the negative of the 
distance between the wall and the second layer. At low 



p, however, there is a rapid ~ 1/p growth and a large Chemical Society, 
slip. 

The slip length is a sensitive function of the properties 
of the wall. When ^4=0, there is large slip even at high 
fluid densities, which is qualitatively consistent with the 
scenario of specular collisions with the wall. It is also 
reminiscent of the large slip found in MD studies of a 
nonwetting droplet on a solid surface |Ie) | However, for 
intermediate values of A there is a complicated crossover 
behavior, as illustrated by the £ - p curve for A = j and 
|. In Ref. Q, a universal relationship between ( and 
shear stress was found in the dense fluid regime on vary- 
ing V w f . Our large p data on £ can similarly be collapsed 
by rescaling but our control of V w t is through the param- 
eter A as opposed to the amplitude of the full potential. 
The latter situation || would lead to a vanishing of the 
channel confinement in the zero amplitude limit. 

The right-hand panel of Figure 4 shows that the variation 
of the viscous slip length ( with p for chain molecules 
is more sensitive to A than for monatomic fluids. For 
A = 1, £ becomes most negative as p tends to zero, while 
for A — j, C is large and positive even at high p (20.6er 
for p = 0.79er~ 3 ). This result is qualitatively consistent 
with the MD results |L6| on hexadecane in the vicinity of 
a solid surface, and with experiments on thin polymeric 
film drainage |l7| ]. 

For monatomic fluids, the fluid velocity in the middle of 
the channel in a gravity driven flow exhibits a maximum 
as a function of p fn]| . In contrast, the bottom panel of 
Figure 2 shows that, for chain molecules and A—l, the 
dependence of v max on p is monotonic even though the 
viscosity at low p is substantially similar for the n = 1 
and n = 10 fluids (data not shown). The large p behav- 
ior of v max is dictated by the viscous effects associated 
with a dense fluid. The increase in v max with p at low 
p for the monatomic case originates in the active recon- 
struction of the second layer. In contrast, for the chain 
molecule fluid, the organization near the walls remains 
largely unchanged (the chain molecules of the first two 
layers are stuck to the wall) so that viscous effects dom- 
inate the behavior of v max even at small densities. 

In summary, the two modes of behavior of £ as p de- 
creases to a small value, rapid growth or approach to a 
constant, are seen to be associated with distinctly differ- 
ent dependences of v max on p. They also appear to be 
linked to the presence or absence of a density buildup 
in the second layer. Strikingly, the apparent complexity 
and non-universality of the BC is due to the microscale 
fluid physics in the vicinity of the solid interface. 
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FIG. 2. The top two panels show the density, pi<J Z , in 
the first (i=l, square symbols) and second (i = 2, circular 
symbols) layers for the indicated values of A. The top 
and middle panels are for the simple and polymeric fluids, 
respectively. The bottom panel shows the velocity of 
the polymeric fluid in the center of the channel vs. p, 
in gravity driven flow with g = COle/mer. Lq refers to a 
channel of standard size and 2Lq corresponds to channels 
with twice the width. In all cases, the statistical errors 
are negligible. 

FIG. 3. Top: Velocity profiles in Couette flow for var- 
ious p, in the narrowest channel with attractive walls; 
the upper and lower walls move with velocities ±0.1<7/ro, 
respectively. The slip length for the high p case is - 
1.7a. Bottom: The temperature profile in the thermal 
slip problem. The open and closed circles correspond to 
p=0.773 and 0.004 c^ 3 respectively. For p=0.193cr" 3 , 
only the central slope is indicated. 

FIG. 4. Q (measured in units of a) vs. p for viscous and 
thermal slip phenomena, for various values of A. The 
horizontal lines labeled I and II indicate the locations of 
the first and second layers, respectively. In all panels, 
triangles correspond to a channel of doubled width, and 
the remaining symbols are for the standard width. In 
the left and right panels, the solid and open symbols 
correspond to £, as determined from Couette and gravity 
driven (acceleration 0.01e/mcr), flows respectively. 



Figure captions 



FIG. 1 Density profiles for the monatomic fluid under- 
going thermal slip for walls with A=l. The tethering 
centers of the molecules comprising the inner wall are 
at the edges of the figure, and the values of density at 
the center of the channel, p, are indicated next to the 
histograms. The peak which is the closest to the wall 
denotes the first layer, and the second corresponds to 
the second maximum plus the elevated bin next to it. 
The inset shows the difference in the densities of the sec- 
ond layers in the vicinities of the cold and hot walls (the 
subscripts 2c and 2h, respectively) as a function of the 
central density. 
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